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Identified charged pion, kaon, and proton spectra are used to explore the system size dependence 
of bulk freeze-out properties in Cu+Cu collisions at s/snn~ 200 and 62.4 GeV. The data are studied 
with hydrodynamically-motivated Blast-wave and statistical model frameworks in order to charac- 
terize the freeze-out properties of the system. The dependence of freeze-out parameters on beam 
energy and collision centrality is discussed. Using the existing results from Au+Au and pp collisions, 
the dependence of freeze-out parameters on the system size is also explored. This multi-dimensional 
systematic study furthers our understanding of the QCD phase diagram revealing the importance 
of the initial geometrical overlap of the colliding ions. The analysis of Cu+Cu collisions, which 
expands the system size dependence studies from Au+Au data with detailed measurements in the 
smaller system, shows that the bulk freeze-out properties of charged particles studied here scale 
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with the total charged particle multiplicity at mid-rapidity, suggesting the relevance of initial state 
effects. 



PACS numbers: 25.75.-q 

I. INTRODUCTION 

An experimental study of relativistic heavy ion col- 
lisions augments our understanding of the QCD phase 
diagram The high energy density reached in such 
collisions at the Relativistic Heavy Ion Collider (RHIC) 
is believed to result in a novel state of hot and dense 
matter with properties strikingly different from that of a 
hadron gas or ordinary nuclear matter Q. 

The bulk properties of particle production are stud- 
ied using identified particle spectra at low momentum. 
Model-dependent interpretations of the measured data 
provide insight into the complex dynamics of the colli- 
sion and further explore the QCD phases through which 
the collision evolves. The dense system formed in the 
early stages of the collision continuously expands and 
cools, until kinetic freeze-out, beyond which the parti- 
cles stream freely into the detector. Through measure- 
ments of species abundances and transverse momentum 
distributions, information about the final stages of the 
collision evolution at chemical and kinetic freeze-out can 
be inferred. 

The relative particle abundances and spectral shapes 
discussed here were tested within the frameworks of 
statistical (chemical equilibrium) H and Blast-wave Q 
models. In the chemical equilibrium model, particle 
abundances relative to the total system volume (assumed 
to be the same for all particle species) are described 
by system temperature at freeze-out, the baryon and 
strangeness chemical potentials, and the strangeness sup- 
pression factor. The Blast-wave model describes spec- 
tral shapes assuming a locally thermalized source with a 
common transverse flow velocity field. The results from 
Au+Au collisions at v /s^=200 and at 62.4 GeV 0- 
have shown that the chemical freeze-out temperature, 
T c h, has little dependence on centrality whereas the ki- 
netic freeze-out temperature, Jkim decreases with in- 
creasing centrality of the collision. Further, the radial 
flow velocity, /3, increases with increasing centrality. The 
observed changes in Tki n and (3 with centrality are con- 
sistent with higher energy and pressure in the initial 
state for more central events. On the other hand, the 
centrality independence of the extracted chemical freeze- 
out temperature indicates that even for different initial 
conditions, collisions always evolve to the same chem- 
ical freeze-out. Moreover, the value for the chemical 
freeze-out temperature in Au+Au is close to the criti- 
cal temperature, predicted by some lattice QCD calcula- 
tions [3] . This suggests that chemical freeze-out coincides 
with hadronization and therefore T c h provides a lower 
limit estimate for the temperature of the prehadronic 
state @. The systematic behavior of the kinetic freeze- 
out properties with charged hadron multiplicity appears 



to follow the same trend for all energies and systems at 
RHIC d, H|. In this paper, the systematic studies of the 
QCD phase diagram from heavy ion collisions are en- 
riched by the addition of new RHIC data from Cu+Cu 
collisions at v / J/v~/v=200 an d 62.4 GeV. 



II. THE STAR EXPERIMENT 

The Cu+Cu data presented here were collected by the 
STAR experiment during the RHIC 2005 run. Cop- 
per nuclei ( 63 Cu) were collided at y / ijv r /v=200 and 
62.4 GeV. Data were recorded with a minimum bias 
trigger obtained from the Beam-Beam Cherenkov coun- 
ters fit)} coupled with information from the Zero-Degree 
Calorimeters [ll|. This trigger is found to be sensitive 
to the top ~85% of the inelastic cross-section. The data 
studied here correspond to the top 60% of the inelastic 
cross-section (minimum bias) where little or no ineffi- 
ciency of the triggering or vertex reconstruction is found. 
These 0-60% minimum bias events, with 24 M and 10 M 
events recorded at 200 and 62.4 GeV, respectively, were 
divided into six centrality bins each corresponding to a 
10% interval of the geometric cross section. Transverse 
mass distributions for charged pions, charged kaons, pro- 
tons and antiprotons, previously reported by STAR for 
Y / J/vtv=200 GeV pv collisions and Au+Au collisions at 
62.4 and 200 GeV [5ft7|, are used for comparison. 

The STAR Time Projection Chamber (TPC) [H 
tracks particle trajectories over a wide range of momen- 
tum at mid-rapidity (\r]\ < 1.8). The particle identifi- 
cation at 1ow-pt uses measurements of truncated mean 
ionization energy loss, (dE/dx), of the charged particles 
traversing the TPC. Particles of different mass show dis- 
tinct patterns in the (dE/dx) dependence, as shown in 
Fig. [TJ left panel. This allows statistical separation of 
pions and kaons in the momentum range 0.25 < px < 
0.80 GeV/c at mid-rapidity (|y| < 0.1), and of pro- 
tons and anti-protons from other species in the range 
0.40 < p T < 1.20 GeV/c. 

The momentum measurement is given by the curvature 
of the particle trajectories as they pass through the 0.5 T 
magnetic field of the STAR detector. To ensure optimal 
dE/dx resolution, only primary tracks, with dca (dis- 
tance of closest approach between the particle trajectory 
and the event vertex) less than 3 cm, and at least 25 out 
of 45 possible fit points are used in this analysis. Particle 
identification at mid-rapidity (|y| < 0.1) is achieved by 
fits to the Z- variable, defined as a logarithm of (dE/dx) 
divided by the theoretically expected value for each par- 
ticle type, given by Bethe-Bloch formula [llj]. This new 
variable is introduced to remove the strong px depen- 
dence at low momenta. Such a. normalized distribution 
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FIG. 1: The left panel shows the truncated mean ionization energy loss ({dE/dx)) in the TPC as a function of transverse 
momentum for positively charged tracks from 200 GeV Cu+Cu collisions. The right panel shows Z(ir), the logarithm of the 
measured (dE/dx) divided by the theoretical expectation for energy loss of charged pions, for 0.40 < pr < 0.55 GeV/c. Also 
shown is an example four-Gaussian fit that is used to extract the raw yields for different species. 
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FIG. 2: Estimated fraction of background protons in the raw 
proton sample as function of transverse momentum, for the 
most central ^/snn~200 GeV Cu+Cu collisions. No strong 
centrality or energy dependence for this correction was ob- 
served for all Cu+Cu data available. 



is created for a given particle and centrality and is di- 
vided into narrow transverse momentum slices (width 
Apr = 50 MeV/c). These momentum projections are 
fit with a combined four-Gaussian function, one for each 
of the particle species of a given charge: n, K, p and 
e. The integral of each Gaussian provides the raw yield 
at each momentum. This procedure is repeated for each 
particle species in order to assign the correct rapidity for 
each track, using the mass of the particle. Thus, fits to 



the auxiliary particles in each distribution (for example 
K, p and e for tt analysis) are used only to estimate the 
contamination when bands overlap. The right panel of 
Fig. Q] shows an example for pion yield extraction for one 
momentum slice. For more details see Ref. Q- 

The raw yields extracted from each of the four- 
Gaussian fits are then corrected for detector acceptance, 
single-track reconstruction efficiencies, and other effects 
as discussed below. To determine the correction factors, 
simulated tracks were embedded into real data on the 
raw signal level and run through the standard recon- 
struction chain. The estimated single-track reconstruc- 
tion efficiency is about 80% for n in Cu+Cu collisions 
and exhibits a small centrality and pr dependence. The 
PT-spectrum has also been corrected for the energy loss 
by multiple scattering beyond that for pions which is 
calculated during reconstruction. This affects the re- 
constructed momentum at low values. The maximum 
value of this mass-dependent correction to the measured 
Pt value for K ± and p(p) was found to be 2% and 3%, 
respectively, for the lowest measured pt bin. An addi- 
tional correction for the background contamination in the 
proton sample is made. The background protons arise 
predominantly from secondary interactions in the beam 
pipe and detector material (knock-out protons). It is es- 
timated from data to be about 40% at pr = 400 McV/ c, 
diminishing to near zero at px = 1 GeV/ c, as shown 
in Fig. [2] To estimate this correction factor we com- 
pare the distribution of proton dca to that of the anti- 
protons (see Ref. Q for more details). In the measured 
dca region, the integral difference between protons and 
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FIG. 3: The top row shows negatively charged pion (leftmost column), kaon (center) and anti-proton (right) spectra from 
Cu+Cu collisions at y / ijviv=200 GeV. Six centrality classes are shown as dark (central 0-10%) to light (50-60%) shades. 
Spectra for 62.4 GeV Cu+Cu are shown on the bottom row. Statistical and systematic errors (which do not exceed 7%) are 
smaller than the symbol size. 



anti-protons (after normalization by the anti-proton to 
proton ratio) is considered to be the background contri- 
bution to the proton yield. Pion yields are additionally 
corrected for feed-down contributions from weakly de- 
caying particles, muon contamination, and background 
pions from detector material. This correction is found 
to decrease from about 15% at 0.3 GeV/c to about 5% 
at 1 GeV/c. The (anti)protons presented in this paper 
are inclusive measurements (not corrected for weak de- 
cays). It has been found in previous studies that the 
analysis cuts used for the low-p-r identified proton stud- 
ies (dca < 3 cm) reject only a negligible fraction of 
daughter protons from the hypcron decays [Bj]. There- 
fore, our sample reflects the total baryon production in 
the collision. Earlier Au+Au studies [5j and preliminary 
Lambda-hyperon spectra from Cu+Cu collisions [l4[ in- 
dicate that the freeze-out spectral shapes are similar for 
As and protons, resulting in similar spectra shapes for 
primary and feed-down protons. The fraction of the 
weak-decay feed-down protons is estimated to be about 

30% m. 

This analysis technique is used to obtain the low- 
Pt particle spectra for all centrality bins at both 200 and 
62.4 GeV center-of-mass energies and for the Cu+Cu and 
Au+Au colliding systems. Additional technical details 
on the analysis and applied corrections can be found in 



Refs. [!, IH with a thorough overview in Ref. Q- 



III. RESULTS 

The transverse momentum spectra are shown in 
Figs. [3] and @] for tt (leftmost column), (center) and 
(anti)protons (right) in Cu+Cu collisions. The top row 
presents the data for y / s]v r iv=200 GeV, whilst data for 
62.4 GeV are shown in the bottom row. The symbol 
shades represent different centrality bins. The particle 
and anti-particle spectral shapes are similar for all species 
in each centrality bin. At both collision energies a mass 
dependence is observed in the slope of the particle spec- 
tra. Due to the large number of events recorded and good 
tracking efficiency the statistical errors are less than 1%. 
The systematic uncertainties are similar to those deter- 
mined in prior analyses of low-p<r spectra in Au+Au col- 
lisions 0- Systematic errors are divided into two classes: 
point-to-point and scale uncertainties. The overall scale 
uncertainty, mostly due to the embedding procedure for 
the single-track reconstruction efficiency, is estimated to 
be 5% for all particle species. Point-to-point uncertain- 
ties are determined for each px bin and particle species. 
For pions and kaons, this error is evaluated to be less than 
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FIG. 4: The top row shows positively charged pion (leftmost column), kaon (center) and proton (right) spectra from Cu+Cu 
collisions at y / 5jv"iv=200 GeV. Six centrality classes are shown as dark (central 0-10%) to light (50-60%) shades. Spectra for 
62.4 GeV Cu+Cu are shown on the bottom row. Statistical and systematic errors (which do not exceed 7%) are smaller than 
the symbol size. 



7% and 13%, respectively. These maximal errors repre- 
sent pt bins where a significant (dE/dx) overlap occurs 
between tt^ , or . For protons and anti-protons, the 
maximum error is 5%. At \ow-px, the proton uncertainty 
is greater than that for anti- protons (4.0% versus 1.3%, 
respectively, at pr = 400-450 MeV/c) owing to the ad- 
ditional uncertainty from the proton's background. The 
uncertainty due to the background decreases rapidly from 
3.7% at p T = 400-450 MeV/c to 1.5% for p T > 1 GeV/c. 

For the anti-particle to particle yield ratios, systematic 
errors are much reduced due to a cancellation of the effi- 
ciency uncertainties and a partial cancellation of extrap- 
olation uncertainties, as described above. A systematic 
uncertainty of 2%, 3%, and 5% is assigned to 7r~/7r + , 
K _ /K + , p/p, respectively. 

We further fit the obtained pt distributions to extract 
system properties at different stages of the collision evo- 
lution. The first fit to the data probes collision properties 
at kinetic freeze-out. Here, a Blast- wave model Q is used 
to simultaneously fit the 7T , K and (anti)proton spec- 
tra at a given centrality. This fit provides a good descrip- 
tion of the spectra shapes, as illustrated in Fig. [5] with 
results from most central 200 GeV Cu+Cu data. The tt^ 
data points for px < 0.5 GeV/c are excluded from the 
Blast-wave fits to reduce the effects of resonance decay 
contributions as done in previous works Including 



this low-pT region in the fit leads to a poorer description 
of proton and kaon shapes, however the resultant modi- 
fication of the extracted parameters remains well within 
their systematic uncertainty. The freeze-out parameters 
obtained from this model are discussed later. Also shown 
in Fig. [5] are Bose-Einstein (cx l/(exp — 1)) fits to the 
7T ± , which provide a slightly better description of these 
data. For evaluation of the systematic uncertainties from 
extrapolation, m^-exponential (cx 1/ exp ) and Boltz- 
mann (cx tut / exp ) fits are also used in the analysis 
(for more details see Ref. @)- 

The particle mean-pr and total particle yields at mid- 
rapidity (|y| < 0.1) are shown in Fig.[S]and Fig. respec- 
tively. The values presented for kaons and (anti)protons 
are determined from the measured spectra points, extrap- 
olated outside the fiducial range using Blast- wave fits dis- 
cussed above. Similarly, a combination of the measured 
data-points and extrapolation from Bosc-Einstcin fits is 
used for the pions. The measured fraction of the total 
yield is found to be 62% for tt ± , 58% for K± and 65% for 
(anti)protons for the most central 200 GeV data; these 
fractions are slightly higher in other centrality bins and at 
lower energy [16| . The systematic uncertainty on dN/dy 
and mean-py, shown in the figures, includes the extrap- 
olation uncertainty evaluated by means of various model 
fits mentioned earlier. Overall, these are estimated to 
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be near 15% of the yields outside the fiducial range for 
pions and kaons and 15%-25% of the extrapolated yields 
for protons and anti-protons, depending on ccntrality. 

We also determine the total charged hadron production 
per unit of pseudo-rapidity, dN^/dr], at mid-rapidity. 
The total particle yield at mid-rapidity for each species, 
obtained by extrapolating the fits to the measured spec- 
trum in the momentum range outside our fiducial cov- 
erage, was corrected for the Jacobian transformation 
dNch/dy—>dNch/dr). The sum of the total charged pion, 
kaon and (anti)proton yields was then corrected for the 
feed-down of weakly decaying neutral strange particles, 
providing the estimate of primordial charged hadron 
yield at mid-rapidity. A complementary method was also 
used, integrating over the charged hadron spectra cor- 
rected for efficiency, feed-down and the Jacobian trans- 
formation, and yielded consistent results. 

The mean-pr of each particle species (ir, K, p) in- 
creases with the number of charged hadrons at mid- 
rapidity dN c ii/dT], as shown in Fig. [BJ Moreover, the 
mean-pr for each particle species appears to scale with 
dN^/drj at mid-rapidity, and to be independent of the 
colliding system and the center-of-mass energy. The par- 
ticle yields show the same systematic scaling features 
with dN c h/dr] as mean-p^ across system and collision 
energy, see Fig. [7] In this logarithmic representation the 
particle yields for each species appear to increase linearly 
with multiplicity, with Cu+Cu matching the Au+Au 
data at similar values of dN^/drj. When shown on a 



linear scale, the integrated yields exhibit a near-linear de- 
pendence with dN c h/dri. The logarithmic scale for both 
axes used here preserves the apparent linear dependencies 
whilst better illustrating the lower multiplicity Cu+Cu 
data. For the detailed features we investigate the rela- 
tive particle production in the following. 

The relative abundances of particles provide an impor- 
tant insight into the chemical properties of the system. 
The relative kaon yield reflects the strangeness produc- 
tion in the collision, whereas proton with respect to pion 
production is dependent on the baryon production and 
transport. Figure [SJa) shows the ratios for the nega- 
tively charged particles, p/n~ and K~/7r - , as a func- 
tion of dNch/dr], which exhibit similar d/V c h/G?77-scaling 
behavior at each collision energy. The slight decrease 
of the values for both ratios seen at the lower collision 
energy of 62.4 GeV is insignificant within experimental 
uncertainties. Figure [8^6) shows the ratios for positively 
charged particles, p/ir + and K + /ir + , which also exhibit 
a diV c h/d?7-scalmg behavior within the same collision en- 
ergy. The beam-energy effect is reversed here as com- 
pared to the ratio of negatively charged particles. Sum- 
ming over the two charges (Fig. Hfc)), the corresponding 
ratios exhibit a common scaling behavior with dN^/dr], 
independent of colliding system and collision energy. The 
energy dependence of the positive and negative parti- 
cle ratios considered separately, points to the effects of 
baryon transport to mid-rapidity, which decreases with 
increasing energy (see also [l6|). 
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Wc further explore the kaon production in Cu+Cu col- 
lisions to gain better insight into production of strange 
quarks. The K~ /ir~ ratio scales with dN^/drj at both 
energies and there is no hint of an additional strangeness 
enhancement of charged kaons in the smaller Cu+Cu 
system compared to the larger Au+Au system. Early 
works from SPS energies reported such additional rela- 
tive strangeness enhancement in the K/ir ratio for smaller 
systems, althou gh no final confirmation of this observa- 
tion is available |l7l . [l8| . The pion and kaon enhancement 
factors arc compared in Fig. [5] This factor is defined 
as the yield per mean number of participating nucleons 
(estimated using a Glauber model), iV part , in heavy-ion 
collisions divided by the respective value in pp collisions. 
A progressive enhancement of kaon production with re- 
spect to pions as a function of collision centrality is evi- 
dent, as shown earlier by the K/ir ratios (Fig. [5{a) and 
Fig. |8]^)). A comparison of these enhancement factors 
between Cu+Cu and Au+Au data is also shown. The en- 
hancement factors for kaons do not show universal scaling 
features with respect to Af part , and are indeed found to 
be higher in Cu+Cu collisions compared to the Au+Au 
system. However, these features do not appear to be 
unique to kaons. A similar trend is observed in Fig. [5] in 
the pion enhancement factors for the two systems. This 
suggests that the additional enhancement, seen in the 
charged kaon yields, is not related to strangeness produc- 
tion, but other physics mechanisms, for example, addi- 
tional entropy production. It should be noted that while 



FIG. 7: Integrated yields at mid-rapidity for pions, kaons 
and anti-protons as a function of the charged particle density 
(dN c h/dri), which is used as a measure of the centrality. For 
comparison Au+Au data are shown as bands. Filled (open) 
points/bands depict data at - v /sjvat=200 GeV (62.4 GeV). 
Error-bars represent statistical and systematic uncertainties 
added in quadrature. 



comparing more spherical central Cu+Cu collisions with 
semi-peripheral Au+Au collisions, the initial conditions 
may not be reflected by N part alone. 

In contrast to pions and kaons, protons show mini- 
mal evolution with centrality and no difference between 
Cu+Cu and Au+Au systems is observed. Figure [5J<i) il- 
lustrates the difference in anti-proton and proton produc- 
tion across energies. As observed in other collision sys- 
tems, the ratio is found to increase and becomes closer to 
unity for higher energy collisions fl9j . The anti-proton to 
proton ratio gives information on the amount of baryon 
transport. In line with the earlier STAR results, our mea- 
surements indicate that while a finite excess of baryons 
over anti-baryons is still present at RHIC energies, p — p 
pair production becomes an important factor. Little or 
no change due to an increase in the system size (central- 
ity) is apparent in the Cu+Cu data at 200 GeV, while 
62.4 GeV data show a decreasing trend with increasing 



9 



or 



0.2 



0.1 



-Cu+Cu 
Au+Au 



0.2 



0.1 



P+P 

-Cu+Cu 
Au+Au 



1.0 
0.8 



P+P 200 GeV 



' pin' 



Kin 



(a) 



200 GeV 
62.4 GeV 
200 GeV 
62.4 GeV 





1 


ffl 




5 


j 




i 

[A A 


I 

A 


A 

a 




i 


i 
















J 



I 

200 GeV 
200 GeV 
62.4 GeV 
200 GeV 
62.4 GeV 



p/n + ' K + /it' + ' 



(b) 




P+P 200 GeV 

n ? UCu+Cu 200 GeV 
0.2 |-ou+uu 62 4 GeV 

200 GeV 
62.4 GeV 



0.1 



(p+p)/jt K/Ji 



(c) 



Au+Au 



P+P 200 GeV 

^Cu+Cu 200 GeV 
hOu+t,u 62 4 QeV 

Au+Au 200 GeV 

AU+AU 62A GeV 



^pTp" 



(d) 




dN ch /dn 
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rapidity. Error-bars represent statistical and systematic un- 
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centrality for this ratio for both Cu+Cu and Au+Au 
data. 



IV. FREEZE-OUT PROPERTIES 

The particle yields and their ratios provide further in- 
formation on the thermal properties of the system at ki- 
netic and chemical freeze-out. 



A. Kinetic Properties 

The completion of all elastic scattering marks the final 
stage of collision evolution and could be interpreted as a 
kinetic freeze-out, where the particle momentum spectra 
are fixed. To quantify this stage, fits are made simultane- 
ously to the spectra of all particle species, but indepen- 
dently for each centrality class (see Fig. [S]for example). 
The fits used here are based on the previously discussed 
Blast- wave model which assumes a radially boosted 
thermal source. These hydrodynamically-motivated fits 
describe the mass dependence of particle spectral shapes 
in terms of the radial flow velocity (/3) , the kinetic freeze- 
out temperature (Tki n ) and the flow velocity profile expo- 
nent (n) at the final freeze-out. The extracted value for 
n is not used to derive any physics interpretation. The 
effects from resonance contributions to the pion spectral 
shape are reduced by excluding the low-p^ data points 
(below 0.5 GcV/c). To enable a comparison with ear- 
lier results on pp and Au+Au collisions 0, Q , the same 
model and the same procedures for the fits are adopted, 
thereby avoiding any possible systematic bias. 

The Blast- wave fit results for the temperature of freeze- 
out are shown in Fig. I10I Tki n and (3 show similar depen- 
dences as a function of dN^/dr] in both Cu+Cu and 
Au+Au collisions, evolving smoothly from the lowest 
to the highest multiplicity, from pp to central Au+Au. 
21in decreases smoothly with centrality implying that 
freeze-out occurs at a lower temperature in more cen- 
tral collisions. The similarity of kinetic freeze-out pa- 
rameters in the events with similar multiplicity from dif- 
ferent colliding species is confirmed by the data alone. 
As noted earlier, the particle mean-pr increases with in- 
creasing dN c h/dr), which is consistent with an increase of 
radial flow with centrality. We note, however, that other 
physics mechanisms, for example, hard and semi-hard 
scatterings, can contribute to higher mean-pr values ob- 
served for kaon and proton spectra [2(|. Direct spectral 
shape comparisons of Cu+Cu and Au+Au events from 
similar multiplicity bins, shown in Fig. [Til show the same 
PT-dependencies between pion spectra from the two sys- 
tems. The same is seen to hold for the respective kaon 
and proton spectra. The middle panel of Fig. |TD] shows 
in addition the chemical freeze-out temperatures for dif- 
ferent colliding systems at different energies. Both the 
chemical freeze-out and the kinetic freeze-out tempera- 
ture show similar scaling features, reflecting the common 
trends in mcan-px and the ratios oip/n and K/-7T, dis- 
cussed earlier. Similarly, on the left panel of Fig. [TU] we 
observe a common dA^h/c^-dependence for the average 
radial flow velocity at kinetic freeze-out. 

A more important observation is that the obtained 
kinetic freeze-out parameters for pions, kaons and 
(anti)protons follow the same trends with dN^/drj, in- 
dependent of collision energy, even though the produc- 
tion cross-sections of the underlying spectra are differ- 
ent. This observation suggests that the kinetic freeze-out 
properties are determined by the initial state. Further- 
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FIG. 9: Enhancement factors for negatively charged pions (left), kaons (center), and anti-protons (right) as function of iVpart in 
y / iiviv=200 Cu+Cu and Au+Au collisions. Error-bars represent statistical and systematic uncertainties on the A+A measure- 
ments added in quadrature. The shaded bands depict model uncertainties on number of participants calculation. The bands 
on the left show uncertainties from the pp measurements that are correlated for all data points. 




FIG. 10: Comparison of kinetic freeze-out properties obtained from fits to Cu+Cu (symbols) and Au+Au (bands) collision data 
at ^/snn~200 (closed symbols/bands) and 62.4 GeV (open). The kinetic freeze-out temperature, Tkin, is shown versus flow 
velocity, /3, and multiplicity in panels (a) and (b) respectively (more central collisions are to the right side of each plot). Panel 
(b) also shows the multiplicity dependence of the chemical freeze-out temperature, T c h, (square symbols). Panel (c) shows the 
multiplicity dependence of the average radial flow velocity. 



more, a modcl-dcpcndent connection between the num- 
bers of produced charged particles and the initial gluon 
density of the colliding system [2l| can be used to deduce 
that the freeze-out properties are most probably deter- 
mined at the initial stages of the collision and are driven 
by the initial energy density. 



B. Chemical Properties 

Chemical freeze-out occurs at the stage of the collision 
when all inelastic interactions cease and the produced 
particle composition in terms of yields is fixed. Valu- 
able information for this collision stage can be obtained 
directly from the experimental results by forming parti- 
cle ratios and comparing them across different collision 
systems and energies. 



The ratios of the different particle yields in Cu+Cu col- 
lisions arc further analyzed within the framework of the 
statistical model 0]. This model describes the chemi- 
cal freeze-out of the colliding system by several fit pa- 
rameters: the temperature at which freeze-out occurs 
(T c h), the cost of producing matter in terms of baryon 
and strangeness chemical potentials (hb, Ms); and an 
additional ad-hoc parameter, known as the strangeness 
suppression factor, (7 S ), to reconcile the lower yield of 
strange hadrons in collisions involving smaller species (for 
example pp and d+ Au) . 

These statistical fits are performed on the relative par- 
ticle abundances from tt^, and p(p) alone. Figure [T2l 
shows an example of the resultant fit to the identified 
hadron ratios from central ^/sNN=200 GeV Cu+Cu col- 
lisions. The lower panel of this figure illustrates the fit 
quality. We note, that the successful description of the 
ratios by the model could not prove the attainment of 
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FIG. 11: Comparison of the spectral shape between Cu+Cu and Au+Au data at y / ijviv=200 GeV. Centrality classes are chosen 
with a similar average charged hadron multiplicity at mid-rapidity. Pion (left), kaon (center) and anti-proton (right) spectra 
are shown for 10-20% central (40-50%) Cu+Cu (symbols) compared to 40-50% mid-peripheral (60-70%) Au+Au (lines). 
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chemical equilibrium, but suggests the statistical nature 
of particle production in these collisions p3 |. The re- 
sults obtained for the freeze-out parameters are shown in 
Figs. ES1GJ and M 

Statistical model fits to a wider variety of hadron yields 
were also attempted using preliminary results for the A, 3 



FIG. 13: Baryon and strangeness chemical potentials, fiu and 
[is, as a function of dN^/dr) for 200 and 62.4 GeV in Cu+Cu 
(symbols) and Au+Au collisions (bands). 



and 4> particles and anti-particles from 200 GeV Cu+Cu 
data from [l4l ]. Including more particles into the model 
fits reduces the systematic uncertainty on the extracted 
parameters and resulted in parameter values consistent 
with those obtained from fits using tt , K ± and p(p) 
alone reported here. In general, the observed systematic 
trends in the freeze-out parameters as a function of the 
collision centrality are preserved 0, 0] • 

Figure Q3] (left panel) shows the evolution of the chem- 
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suppression factor, 7s, as a function of dN c h/drj for 200 and 62.4 GeV in Cu+Cu (symbols) and Au+Au collisions (bands). 



ical freeze-out temperature versus baryon chemical po- 
tential in central heavy-ion collisions from the very low 
energy SIS data through AGS and SPS to RHIC (STAR 
data points only) . The overall evolution of T c h can be re- 
produced by the phenomenological model fit [23| applied 
here to all the data points shown (dashed line). As the 
collision energy increases, the temperature at freeze-out 
is found to increase up to SPS energies. This is followed 
by a plateau at RHIC energies at a value close to that 
of the hadronization temperature expected from lattice 
QCD calculations. At RHIC, for all systems and center- 
of-mass energies, T c h appears to be universal, as shown 
in Fig. [10] (middle panel) . 

The value of the baryon-chemical potential at a given 
center-of-mass energy is found to be slightly higher for 
the larger system, with Au+Au and Cu+Cu measure- 
ments showing common trends with charged hadron mul- 
tiplicity (Fig. [1_2]). We note that, presented in the same 
figure, values of strangeness chemical potential are close 
to zero with no obvious systematic trends for all ener- 
gies and colliding systems studied at RHIC. Within a 
given system, (1b reflects the decrease in net-baryon den- 
sity with increasing collision energy from ^/sjv./v=62.4 to 
200 GeV. This behavior can be observed directly from 
the particle ratios, where p/p increases as a function of 
energy (Fig. [5]). For the most central Cu+Cu events we 
measure p/p = 0.80 ± 0.04 at 200 GeV, and 0.55 ± 0.03 
at 62.4 GeV. The lack of centrality dependence in the 
baryon to meson ratios in Cu+Cu and Au+Au data, 
points to similar freeze-out temperatures for the studied 
systems. The constant values of T c h at RHIC energies 



for collisions with different initial conditions, energy and 
net-baryon density, points to a common hadronization 
temperature of the systems. 

Another parameter extracted from the fit, which is re- 
lated to strangeness production, is the strangeness sup- 
pression factor, 7s, shown versus dN^/drj in Fig. 1141 
The suppression of strange hadron yields is observed in 
smaller systems, such as pp and peripheral collisions. 
Within statistical models this can be explained by a re- 
duced production volume 24j. At low beam energies, 
where equilibration of s quarks with respect to u and d is 
not expected, the suppression is also seen. We find that 
within the systematic errors on the fit parameters the 
strangeness suppression factor in Cu+Cu is consistent 
with that for Au+Au for the same number of charged 
particles, dN^/dr]. As only charged kaon yields were in- 
cluded in the fit, this observation is directly related to 
an absence of any additional enhancement in K /ir at the 
same dN^/dr] in the smaller Cu+Cu system with respect 
to the larger Au+Au system as discussed previously 

The 7 S parameter shows a similar increase with cen- 
trality for both systems and energies. The value of 7 S 
approaching unity for the central Au+Au collisions in 
the context of thermal model would imply that the pro- 
duced strangeness is close to equilibrium. 



V. SUMMARY 

We have presented measurements of identified charged 
hadron spectra in Cu+Cu collisions for two center-of- 
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mass energies, 200 and 62.4 GeV. These new results 
of 7r , and p(p) have further enriched the variety 
of 1ow-pt spectra at RHIC. The data have been stud- 
ied within the statistical hadronization and Blast-wave 
model frameworks in order to characterize the properties 
of the final hadronic state of the colliding system as a 
function of system size, collision energy and centrality. 

These multidimensional systematic studies reveal re- 
markable similarities between the different colliding sys- 
tems. No additional enhancement of kaon yields with 
respect to pions is observed for the smaller Cu+Cu sys- 
tem compared to Au+Au. The obtained particle ratios, 
mean-pT and the freeze-out parameters, including the 
strangeness suppression factor, 7 S , are found to exhibit 
a smooth evolution with dN c h/dr], and similar proper- 
ties at the same number of produced charged hadrons 
are observed for all collision systems and center-of-mass 
energies. The bulk properties studied have a strong cor- 
respondence with the total particle yield. Within thermal 
models this reflects a relation between the energy per par- 
ticle at freeze-out and the entropy derived from particle 
yields, which reflects the initial state properties for adia- 
batic expansion. The baryon chemical potential could 
in addition be influenced by the initial valence quark 



distribution and by baryon transport during expansion, 
leading to a more complicated dependence. The scal- 
ing features of freeze-out properties are not presented at 
the same A part for lighter and heavier ions as scaling is 
badly broken when data measured at different energies 
are compared. This suggests that N part does not reflect 
the initial state of the system accurately. 
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